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We discuss the nonlinear current of an interacting quantum dot coupled to normal and superconducting reser¬ 
voirs with applied voltage and temperature differences. Due to the particle-hole symmetry introduced by the 
superconducting lead, the pure (subgap) thermoelectric response vanishes. However, we show that the Andreev 
bound states shift as the thermal gradient increases. As a consequence, the I—V characteristic can be tuned 
with a temperature bias if the system is simultaneously voltage biased. This is a cross effect that occurs beyond 
linear response only. Furthermore, we emphasize the role of quasiparticle tunneling processes in the generation 
of high thermopower sensitivities. 

PACS numbers: 74.25.fg, 74.45,+c, 74.78.Na, 73.23.-b 


I. INTRODUCTION 

When a superconductor (S) material is attached to a nor¬ 
mal (N) conductor the transfer of charges across the NS in¬ 
terface is dominated by Andreev processes Oh in which elec¬ 
trons (holes) are retroreflected as holes (electrons) adding (de¬ 
stroying) Cooper pairs into the ground state. Leakage of su¬ 
perconducting pairing correlations into the normal conductor 
has a profound impact when sandwiched quantum dots (QDs) 
are considered. Hence, multiple Andreev reflections lead to 
the formation of Andreev bound states observed in N-QD-S 
tunnel experiments when the normal contact acts as a probe 
terminal B-'S]- These hybrid setups are excellent testbeds to 
examine the interplay between Coulomb interaction and prox¬ 
imity effects. In particular, N-QD-S systems exhibit remark¬ 
able conductance changes in the cotunneling regime @] orig¬ 
inated from the occurrence of zero-bias anomalies, which at 
much lower temperatures arise from the competition between 
Kondo physics 0] and the Yu-Shiba-Rusinov states (usual 
Andreev bound states in magnetically active platforms). In¬ 
deed, the Yu-Shiba-Rusinov states are viewed as precursors 
of Majorana quasiparticles in hybrid systems with large spin- 
orbit interaction and magnetic fields iMlSll- Additionally, N- 
QD-S structures have been proposed as suitable sources of 
solid-state qubits [ 10). 

Previously cited works cope solely with the electric re¬ 
sponse of hybrid systems. In contrast, their thermoelectrical 
response is much less understood. In QDs coupled to nor¬ 
mal electrodes, the thermoelectric voltage V* generated in re¬ 
sponse to a small thermal gradient AT is greatly amplified 
across each dot resonance [1 fi ll 2I1 . leading to large values of 
the Seebeck coefficient S = —V t h/AT ill 3:1. Superconductiv¬ 
ity can significantly alter these thermoelectric properties El 
and provide additional information, as demonstrated with An¬ 
dreev interferometers El El. Importantly, in hybrid nanos¬ 
tructures large values of S have been envisaged by breaking 
the particle-hole symmetry ElS- Here, we focus on N- 
QD-S setups that indeed preserve such symmetry and thereby 
avoid the generation of electrical currents by thermal gradients 
alone. Remarkably enough, we find that the nonlinear trans¬ 
port regime do exhibit a large cross thermoelectric coupling 
signal when both voltage and temperature shifts are present. 



FIG. 1. (Color online) Top: A quantum dot device tunnel coupled to 
hot normal (N) and cold superconducting (S, gap 2A) leads whose 
distribution functions are sketched. Andreev processes where in¬ 
jected charge forms a Cooper pair in the S lead are indicated for 
(a) incident electrons and (b) incident holes. For voltage bias V = 0 
both processes counterbalance and the net Andreev current vanishes 
even in the presence of a thermal gradient 9. Therefore, a cross 
coupling (finite V and 9) is needed to generate thermoelectric trans¬ 
port. Bottom: Andreev transmission for different thermal bias 9. 
Tilted arrows indicate Andreev level shifts of the bound states indi¬ 
cated in (a,b) when the thermal bias enhances from ksd = 0.1 A to 
ks9 = 0.5A. Parameters: £d = 0.15A, Tl = 0.1A, Tr = 0.5A, 
ksT = 0.1A, V = 0. 

Below, we discuss the details. 

The nonlinear regime of transport have been hitherto poorly 
explored in interacting hybrid systems (for an exception, see 
Ref.[2db. The nonlinear thermoelectric transport is unique be¬ 
cause the conductor nonequilibrium potential responds in a 
nontrivial way to large shifts in the applied thermal ( 9 ) and 
voltage ( V ) differences [ 21 ~ [29ll . Consider the charge current 
I for a generic conductor, 

I = G 0 V + GxV 2 + L o 6 + L 1 0 2 + M 1 V0 1 (1) 

expanded up to leading order in the rectification terms. For 
the present setup [see top panel in Fig. |T). all pure thermo- 
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electrical coefficients [Lq and L\ in Eq. <|T]) 1 vanish in the sub¬ 
gap region due to symmetry between the processes indicated 
in Fig. ©a) and Fig. [U b), canceling the net current. Impor¬ 
tantly, whereas / is not sensitive to a thermal gradient alone, 
it does react to a simultaneous change in voltage and temper¬ 
ature shifts because for finite V and 9 the process sketched in 
Fig. ©b) is Pauli blocked for positive voltages. This unique 
cross effect is quantified to leading order by the coefficient 
M\ = d 2 I/dVdO | eq (eq means {V = 9 = 0}). This mixed 
term has thus far received little attention, because its effect is 
masked in many systems by the diagonal responses Gi and 
L\. We propose that hybrid systems are highly appropriate 
experimental setups to test this cross response as they exhibit 
a clear signal of M\. Our results show that the I—V charac¬ 
teristic can be thermally tuned due to M\. Furthermore, the 
cross response is intimately linked with the fact that Andreev 
bound states shift their positions [Fig. H]c)] under the influ¬ 
ence of a temperature shift 9, a feature easily accessible with 
tunnel spectroscopic techniques. 

The paper is organized as follows. In Sec. QI| we explain 
the gauge-invariant nonlinear thermoelectric transport theory 
to describe our N-QD-S junctions. Section [HI] is devoted to 
displaying our main results based on this theory. Importantly, 
the cross thermoelectric coupling term Mi in Eq. © plays 
a significant role to the Andreev current I a as illustrated in 
Fig. [3] where I a can increase or decrease its amplitude ac¬ 
cording to the quantum dot energy level. Furthermore, we 
show that high thermopower can be created when combined 
with the quasiparticle contribution uniquely in the nonlinear 
transport regime. Finally, we conclude our findings in Sec.lIVI 


II. HYBRID CONDUCTORS: NONLINEAR 
THERMOELECTRIC TRANSPORT THEORY 


We consider a N-QD-S system with a heated normal lead 
as in Fig. ©a). The model Hamiltonian can be written as 

H = 'Hl+'H r +H d +'Ht, (2) 

where Hl = Yka £ Rkc\ ka CLka describes the electron with 
momentum k and spin cr in the left (L) normal lead while 

Wr = Spo- £ RpC l Rpa c Rp° + £ P [^ c fi,-pt c Rp4. H- c -\ d e " 
scribes electrons with momentum p in the right ( R) supercon¬ 
ducting contact. The second term in Hr describes the Cooper- 
pair formation with energy cost A. The dot Hamiltonian is 
Hd = Y a ( £ d ~ eU)dl r d a , with e d the dot energy level renor¬ 
malized by the internal potential U that accounts for electron- 
electron repulsion. This interaction term will be considered 
at the mean-field level, which is a good ap prox imation for 
metallic dots with good screening properties IBOIl . In Eq. ©, 

Ut = Y,ka t Lc\ k!7 dv + Y.pa t R e * eVRtc R P v d <? + Hc ■ de¬ 
picts the tunneling processes between the dot and each lead 
with amplitudes ip, and t R (Vr is the voltage in lead S). 

We use a gauge-invariant current-conserving theory applied 
to nonlinear thermoelectric transport. First, the electric cur¬ 
rent can be evaluated from the time evolution of total electron 


number in the left lead. Nr = Yka c Lka c Lka, through 

I = -e{N L {t))=-(ie/h)([U,N L ]). (3) 

Within the nonequilibrium Keldysh-Green function formal¬ 
ism mm, one finds that / = I a + Iq is a sum of two 
terms: the Andreev current I a and the quasiparticle contribu¬ 
tion Iq, 


I A = ^JdeT A (e) [f L (e - eV) - f L (e + eV )], (4) 

I Q = Y J de T Q( £ ) [/i( £ “ eV ) - /fl( £ )] ■ (5) 


Here, the Fermi-Dirac distribution function is given by 
fa=L, R {£ ± eV) = {1 + exp[(e ± eV - E F )/k B T a ]}~ \ 
where the electrode temperature T a = T + 9 a is given in 
terms of the base temperature T and the shift 9 a , the volt¬ 
age bias reads V = Vl — Vr and the Fermi level is taken 
as the reference energy (Ep = 0). The Fermi function f r 
for the superconductor is that of the zero gap state, with the 
gap property included in the S density of states. Quite gener¬ 
ally, the transmissions Ta and Tq are functions of U , which 
depends itself on the applied voltage and thermal bias. 

Both Ta and Tq are expressed in terms of the dot retarded 
Green’s functions G T i:j (e) ( i, j = 1,2) in the Nambu space mb 


GS(e) 


G r n(s) Gi 2 (e) \ 

G r 2 i(e) G r 22 (e) J 


( 6 ) 


whose matrix elements are Fourier transforms of the elec¬ 
tronic part G\i(t,t') = —i@(t — t')({df(t), d|(f')}), the 

hole part G 22 (t,t') = — *0(f — t'){{d^(t), d±(t')}}, and 
those parts that connect both electron and hole dynamics 
G\ 2 (t,f) = -iQ(t - t'XHW.W}). G r 21 (t,t') = 
—i<d(t — i')({dj(i), dj(t')}). Here, the subindex 1 refers to 
the electron sector and 2 to the hole part. The Green’s func¬ 
tions relevant to our work explicitly read G r n (e) = [e - £,/ + 

euLV + czl9+ ^ + 4^^3j-^ r (e)] _1 and 

G\ 2 (s) = G r n (e) 2 f £ ^ A ., A r (e), where A r (e) = [e + e d - 

culV - czlO + ^ ^ v'J ! 1 a ^ ]" 1 and Ul and Zl be 

specified below. For the subgap region |e| < A, one should 
make the substitution \/e 1 — A 2 —> i\JA 2 — e 2 . The contri¬ 
bution of Andreev transmission and that of quasiparticle tun¬ 
neling are respectively given by (32tl 


T A (e)=T 2 L \G r 12 (e)\ 2 , (7) 

~ 7A 

Tq( £ ) = r L r fl (|G^| 2 + |Gy 2 - — Re[G^(G^)*]), 

( 8 ) 


with Tl = 2p\t L \ 2 Y k d ( £ ~ £ Lk), Tr = 2Tr\t R \ 2 Y p H £ - 
e Rp ), and f R = r fl 0(|e| - A)|£-|/\/e 2 - A 2 . Note that T A 
describes the process where an electron (hole) incoming from 
the left lead is reflected as a hole (electron) backward into the 
same lead by producing (destroying) a Cooper pair in the S 
contact. Strictly at zero temperature and in the subgap regime 
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\eV\ < A, the Andreev current is the only contribution to the 
total current. However, quasiparticle poisoning ( Iq ) must be 
also taken into account when temperature or voltage are suf¬ 
ficiently large. The quasiparticle transmission Tq comprises 
the conventional tunnel processes and those in which electron 
or holes in the normal part become quasiparticle excitations in 
the superconducting reservoir, either by keeping the Cooper 
pair number invariant or by creating/destroying pairs [|3 1LI321. 

Substituting Eqs. © and © in Eqs. © and ©, we can 
express the total current I in terms of GC (e). Due to the 
presence of interactions, the Green function depends explic¬ 
itly on the nonequilibrium screening potential U which differs 
for each thermoelectric configuration {V a , 0 Q } a Now we 
discuss how to determine U in a system containing pairing 
correlations. For a weakly nonequilibrium state, the poten¬ 
tial can be expanded SU = [u a V a + z a 9 a ] up to first 
order in the shifts V a and 9 a . Here, SU = U — U eq mea¬ 
sures deviations of the internal potential from equilibrium and 
u a = {dU/dV a ) e q and z a = (dU/d9 a ) eq are the character¬ 
istic potentials describing the system response to the shifts. 
Without loss of generality, we henceforth consider Vl = V, 
Tl = T + 9 (V R = 0, Tr = T). The nonequilibrium poten¬ 
tial satisfies the capacitance equation |33ll 5p = CSU, where 
C is the dot capacitive coupling and the excess charge den¬ 
sity is given by the dot distribution (lesser) Green function 
Sp = p- p eq = if de[G qi {e) - Gf leq (e)]: 




FIG. 2. (Color online) Characteristic potentials: (a) ul and (b) zl 
versus dot level £d at base temperature fcsT = 0.1A with C = 0. In 
(a), we find ul — — 1 for IT 3> Fr (N-dominant) while ul shows 
rather drastic change for IT <Tji (S-dominant). In (b), zl behaves 
qualitatively the same for the two cases. 

Green function: 

= ( |G ii |2 - |C h |2 )«,. (10) 

ZL = era) - 8 -b(i°ni 2 +ig;T)„, 

( 11 ) 

where the integrands are evaluated at equilibrium. Note that 
we have only considered the response to the left lead since 
the superconductor is assumed to be at equilibrium and hence 
urVr + zr9r = 0 irrespective of ur and zr. 


GftOO = ^ |G^| 2 /l(£ - eV) + \G[ 2 \ 2 f L (e + eV) 


r 
27r 


27r . 

/*(e)(|Gn| 2 + |tf 


2A 

121 -^RelG^Gy*]). (9) 


Notice that G[j (e) is dominated by Andreev events when 
energy lies within the subgap region [identified by the term 
/i(e + eV)\Gi 2 {s)\ 2 in Eq. ©] whereas for |e| > A it is 
dominated by (i) the normal electronic dot distribution func¬ 
tion [IG^OOI 2 weighted by left-to-right averaged nonequi¬ 
librium distribution function (Tl/l + r r/r)] and (ii) the 
quasiparticle contribution [last terms in Eq. ©]. Importantly, 
the density is expressed as the sum of injected and screening 
charges, Sp = p m j + p SCI . More explicitly, to leading order in 
V and 9 one has Sp = + D a 9 a ) — U5U where the 

first term in brackets relates to the charge injected from lead a 
with voltage V a and thermal driving 9 a and the last term de¬ 
notes screening effects in terms of the generalized Lindhard 
function n. We can thus define the charge injectivity I) n = 
(dp/dV a ) eq , the entropic injectivity D a = (dp/d9 a ) eq , and 
the Lindhard function n = —( Sp/SU) eq . These quantities, in 
general, contain particle and hole portions (explicit formulae 
are provided in Appendix!© f° r which screening is calculated 
in the presence of superconductivity. 

Solving the capacitance equation, we find SU = 
^2 a {DaVa + D a 9 a )/(C + n), hence we immediately have 
analytic expressions for the characteristic potentials, u a = 
D a /(C + n) and z a = D a j(C + n), in terms of the dot 


III. HYBRID CONDUCTORS: RESULTS 


We next discuss our main findings for the nonlinear ther¬ 
moelectric transport of an interacting hybrid setup. We first 
focus on the internal potential changes SU = urV + zl9 
when the hybrid system is electrically and thermally biased. 
These potentials are shown in Fig. [2] which displays the so¬ 
lutions of Eqs. (flOl) and CD as a function of the dot level 
position for two coupling limits, the normal (T^ > Tfl) and 
superconducting (T^ <C \ r) dominant cases. We consider 
the strongly interacting regime (in the mean-field language) 
by setting C = 0. In Fig. ©a), we find ul — —1 in the N 
dominant case since the applied voltage shifts the dot level 
as £d —> £d ~ ULeV ~ £d + eV due to charge neutrality. 
In contrast, when the superconductor is more strongly cou¬ 
pled to the dot, ul behaves very differently. The sign changes 
(ul > 0) when particle-hole conversion process dominates 
over the ordinary electron tunneling, as expected from the 
term (IG^jJ 2 — |G^ 2 | 2 )eq in Eq. (fl~0t . On the other hand, when 
the dot is thermally driven as shown in Fig.©h) zl shows sim¬ 
ilar behaviors for both coupling limits. The effect of thermal 
driving zl appears only away from the particle-hole symmetry 
point when £d = Ep (recall that Ep = 0), as a consequence 
of a vanishing entropic injectivity. Hence, zl is an odd func¬ 
tion of £d while ml is even. Interestingly, the Andreev trans¬ 
mission Ta is a function of both energy e and potential U and 
thus depends on voltage and 9 via the characteristic potentials. 
This result is illustrated in Fig. |T]'c) for the case of thermal 
driving. Our finding thus suggests an extended controllability 
of N-QD-S junctions using thermoelectric configurations. 
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FIG. 3. (Color online) Andreev current versus voltage for several 
values of the temperature difference ks9 with Fl = 0.1A, Fr = 
0.5A at (a) Ed = 0.2A and (b) £d = 0.7A. In (a), applied thermal 
bias ksd reduces \Ia\ for a given voltage since at Ed = 0.2A, M\ 
in Eq. 0 is negative as shown in the inset of (b). On the contrary, at 
Ed = 0.7A as displayed in (b), 9 increases the amplitude of I a due to 
the positive contribution from Mi. Inset in (a) compares Eq. 0 with 
the exact expression Eq. 0 at zero thermal bias. We take fesT = 
0.1 A and C = 0. 


At low temperature with a small voltage bias \eV\ < A, 
the Andreev process largely contributes to the total current. 
Then, an expansion of Eq. 0 as in Eq. 0 gives the linear 
conductance 


G ° = X/ d£ (-^/) T Ae q , 

the leading-order rectification term 



( 12 ) 


(13) 


and the cross thermoelectric coupling 


Mi = 


4e 2 

h 


ZL- 


dT 'a £ — Ep dT\ 


dU 


T 


de 


■ (14) 


J eq 


As we anticipated above, the thermoelectric response for the 
subgap transport regime (8Ia/ d0)v =o vanishes to any or¬ 
der in 9 , as a result of the Andreev current expression [see 
Eq. 0], Then, L 0 = Li = 0 and Eq. 0 simply becomes 
I = G 0 V + GiV 2 + Mi 1/(9 ( I Q ~ 0 for \eV\ < A). We 
therefore predict that the leading-order thermal response of a 
N-QD-S system is completely determined by the cross ther¬ 
moelectric coefficient Mi. By turning on the electrical bias 
(provided that Ed ^ 0) electron-hole symmetry is lifted and 
then a finite response to a thermal bias is expected. This fea¬ 
ture is unique to Andreev processes since normal tunneling, 
quite generally, gives nonzero Lq and L\ & 

Our finding is now illustrated in Fig. [3] The effect of in¬ 
teractions is clearly visible as I a — V is rectified for nonzero 
values of £d ^ 0: (a) Ed = 0.2A and (b) £d = OTA. Since 
the expansion in Eq. 0 is only valid within the bias range 
|e1//A| <C |G 0 /Gi| and \ks9/A\ <C |G 0 /Mi|, we also 



FIG. 4. (Color online) (a) Quasiparticle and (b) total current versus 
voltage for several ks9 at Ed = 0.3A for Fl = 0.1A, Fr = 0.5A. 
Thermal agitation makes Iq highly contributing to / even for a small 
bias range \eV\, \ks9\ <C A. For a finite ks9 ^ 0 nontrivial solu¬ 
tions of Iq = 0 and 1 = 0 appear (to facilitate the localization of 
Vth from both Iq — V and I — V curves we have drawn a horizontal 
black line at zero y-axis value). We take fcsT = 0.1A. 


compare in the inset of Fig. 0a) with the exact expression 
in Eq. 0. Importantly, the I a- V curves in Fig.0(we choose 
the S dominant case since it offers the clearest signal) can be 
manipulated with the heating gradient 9. The precise behavior 
of I a — 1/ with 9 depends on the dot gate position by either in¬ 
creasing [Fig. 0b)] or decreasing [Fig. 0a)] the amplitude of 
I A with increasing 9. To leading order in V and 9 , this effect 
is explained by the cross coupling coefficient Mi. A salient 
feature of Mi is its dependence of interactions through the 
characteristic potential zl in competition with the noninter¬ 
acting contribution given by the second term in the integrand 
of Eq. ®. The inset of Fig. 0b) shows Mi as a function of 
Ed in the S dominant case, showing that the sign of Mi can be 
tuned by a dot gate potential. This sign determines the lower¬ 
ing or raising of I a — V with 9. We also note that although 
the thermal bias highly affects I a, the current never vanishes 
except for the trivial point V = 0. This differs from the N- 
QD-N case lfl2l [26ll . Hence, a thermovoltage in a N-QD-S 
system can be created via Iq only, an important result which 
we discuss below. 

Unlike the Andreev current discussed above, the expansion 
of the quasiparticle current in Eq. 0 includes all linear and 
lowest order nonlinear coefficients, Iq = GqV + G®V 2 + 
Lq 9 + L®9 2 + M®V9. For low T, Iq is a small contribution 
when \eV\ < A. However, as the thermal bias 9 increases 
Iq starts to contribute nontrivially to the total current I even 
for a small voltage bias \eV\ < A, see Fig. 0a). Importantly, 
we now find finite values of the thermovoltage 14, determined 
from the condition /q(V4) = 0. This is seen in Fig. 0a) as Iq 
left shifts with growing 9, mainly due to nonzero Lq (which 
occurs for Ed / 0). Furthermore, the thermovoltage sign can 
be controlled by the gate potential polarity, i.e., positive (neg¬ 
ative) values of Ed generates negative (positive) 14,. Even if 
the Andreev current cannot exhibit the Seebeck effect solely 
by itself (because no thermovoltage is generated) I a shifts the 
value of 14, when the total current I = Ia + Iq is consid¬ 
ered, see Fig. 0b). Indeed, I a plays a significant role in the 
suppression of thermopower at low thermal biases. 

Figure 0 a) shows the created thermo voltage 14 from the 
total current I for the S dominant case. We find that V4 is 
vanishingly small until a certain amount of 9 is applied, after 
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FIG. 5. (Color online) (a) Thermovoltage versus fcsff/A for several 
Ed- (b) Differential Seebeck coefficient versus base temperature at 
Ed = k B 0 = 0.3A, with T l = 0.1A and r_R = 0.5A. As shown in 
the inset of (a), the linear regime is very narrow suggesting that the 
observed thermoelectric effect is intrinsically nonlinear. Inset of (b) 
displays the thermopower versus the dot level for k B T = 0.1A and 
k B d = 0.3A. 


which Vih boosts. We attribute the suppression of V* around 
/;: n0 < 0.2A to the Andreev current, which lacks the coeffi¬ 
cient of 6 n to any order n as discussed above. We check this 
in the inset of Fig. 0 a), which compares the linear response 
thermo voltage evaluated from I = (Go + G^)V* + L®9 = 0 
with the exact expression from the sum of Eqs. (0) and (0. 
Clearly, the linear approximation quickly fails, implying that 
the thermopower in our system is inherently nonlinear. In 
Fig. 0b), we show the nonlinear Seebeck coefficient S = 
— (dV/dO)j- o at £d = ksO = 0.3A (a moderate value of 
the thermal gradient). We observe a high degree of tunability 
for S' as a function of the background temperature. Another 
tuning parameter is the dot level. The inset of Fig. 0b) shows 
that S sharply increases as Ed is detuned from the particle-hole 
symmetric point. Importantly, the observed thermoelectric ef¬ 
fects can appear with all relevant energy scales well below A. 


IV. FINAL REMARKS 

The combined influence of applied voltages and tempera¬ 
ture biases is a fundamental aspect of electric transport. We 
have here discussed interaction-driven thermoelectric effects 
appearing uniquely in the nonlinear transport regime of a 


IP = - 


IP = - 


Sp p 

m 

Sp h 

w 


/£'-<■> 

/ 2p (£) 


r L 


*|GIi(g)l 

5U 


^|gi 2 (g)| 

su 


For a normal conductor (A = 0) at zero temperature (T = 
0) and close to resonance (e,/ ~ Ep), we obtain the Breit- 


normal-quantum-dot-superconducting system. Andreev pro¬ 
cesses cancel linear Seebeck effects and a nonlinear treatment 
of thermopower is thus called for. We have demonstrated that 
a mixed thermoelectric response determines the thermal driv¬ 
ing of Andeev currents. This result will be robust in the pres¬ 
ence of strong interactions (Coulomb blockade) since even 
when the charging energy is fixed the cross thermoelectric 
coefficient will be nonzero if the transmission is energy de¬ 
pendent, as in most quantum-dot setups. In addition, we have 
found that high thermovoltages can be generated due to quasi¬ 
particle tunneling for moderate thermal gradients, an effect 
which can have crucial importance for superconducting-based 
thermometry and cooling applications 13411 . 
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Appendix A: Injectivities and Lindhard function 


Using the definitions displayed in the main article, we 
determine the charge and entropic injectivities Dl = 
(■ dp/dVi ,)^, Dp = (dp/d9p) eq , and the Lindhard func¬ 
tion n = —(5p/5U) eq . We separate the density function 
(p = p p + p h ) into particle ( p v ) and hole ( p h ) sectors accord¬ 
ing to the Nambu space matrix elements G^ 1 (e) and G\ 2 (e)- 
We find (Dp = D p L +D h L , D L = D p L +D h L , andn = nP+n' 1 ) 


D p - d -f- 

Ul ~ dV 

D h = d JL 

L dV 

d p -2eL 
Ul ~ oe 

nk = d _l_ 
L de 


= ~ erL J ^(-4/)|G^ eq (5)| 2 , 

= eTL j ^(-^) | G Wq( £ )| 2 - 

= -Tl /?^(-4/)|g W£ )| 2 , 


eq 


= -r t 


2t r T 

de e — Ep 
27r T 


(-0 e /)|G; 2ieq OO 


; ^|G! 2 (e) 



eq 


A _S_ 

\T\~8U 


G r i2 [G 


nj 



eq 


Wigner-like approximation for the Lindhard function, n = 

(e/7T)(r/2)[( £d - Epf + (r/2) 2 ]- 1 with r = r L + r R . 
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